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The ground state of three quarks * 
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We measure the static three-quark potential in 5(7(3) lattice gauge theory with improved accuracy, by using 
all available technical refinements, including Liischer-Weisz exponential variance reduction. Together with insight 
gained from 3-state Potts model simulations, our results allow us to sort out the merits of the A- and Y-ansfitze. 
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1. Introduction 

The static qq potential has been studied exten- 
sively on the lattice. In the quenched theory, the 
elementary ansatz V qq (r) = Vq — f + cr qq r turns 
out to be remarkably accurate. The next sim- 
plest system is one of 3 static quarks, from which 
one can gain phenomenological insight about the 
forces inside a baryon. Moreover, the 3 quarks 
must be connected by 3 glue strings to form a 
gauge-covariant object. These strings meet at a 
"gluon junction" , which has been conjectured to 
be a non-perturbative excitation of the QCD vac- 
uum jjj , and might play an important role at the 
hadronization stage in heavy-ion collisions. 

We extract the potential between 3 static 
quarks in SU (3) gauge theory from the exponen- 
tial decay with T of the expectation value of the 
baryonic Wilson loop Fig. 1. We aim at improv- 
ing the accuracy of our earlier work Q through 
several technical refinements. This new accuracy, 
augmented by insight from Potts model simula- 
tions, allows us to reach conclusions about the 
merits of the A-ansatz j3| and the Y-ansatz [Q . 

At short distances, the potential is described 
by perturbation theory, which gives it the form 
(constant + Coulomb), (i) The constant term 
is caused by UV divergences arising from the 
perimeter of the loop. For large T, the bary- 
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Figure 1. Baryonic Wilson loop with junctions at 
x and y:W= ±e abc e a , ¥c ,rfr¥ Tf '. 

lar, mesonic loop, so that the constant term in 
the potential is multiplied by | (for the same lat- 
tice spacing a), (ii) The exchange of one gluon 
between two of the quarks gives \ the one-gluon 
exchange term between a quark and antiquark at 
the same locations. Together, (i) and (ii) imply 



V qqq (n,r 2 ,f 3 ) 



C(r 1 ,r 2 ,r 3 )(l) 
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At large distances, the A-ansatz predicts that the 
potential grows linearly with the perimeter La of 
the quark triangle: V qqq oc a qq ^, so that Eq.(Q) 
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still holds. It is derived from a model of confine- 
ment by center vortices using a beautiful topolog- 
ical argument ||. The F-ansatz predicts instead 
V qqq oc UqqLy, where Ly is the minimal length 
of the 3 flux tubes necessary to join the 3 quarks 
at the so-called Steiner point. It is derived from 
strong coupling arguments Q), and is consistent 
with the dual superconductivity confinement sce- 
nario ||. Since Ly > -^f- for all 3-quark geome- 
tries, the V-ansatz predicts a steeper potential 

Vq Y qq(ri,r 2 ,r 3 ) = V q %(fi,f2,f 3 )+aqq(LY—^)(2) 

with V q ^ q as per Eq.(Q). Both ansatze are con- 
strained to reproduce the diquark limit fj — > ffc, 
V qq q(fi,rj,fk) — ► V q q{?i,Tj) exactly, and there- 
fore contain no free parameter once V q q is given. 
In this respect we differ from the analogous lat- 
tice study of |(5|, where a q q and a qqq are fitted 
separately and therefore not strictly equal. 

2. Technical refinements 

Because the difference between the A- and the 
V-ansatze is very small (I < L ^ 2 < -^=), high 
accuracy in the determination of V qqq is manda- 
tory. The main difficulty at large quark sepa- 
ration is the contribution of excited qqq states. 
Besides smearing the spatial links as in Q, we 
use three additional techniques to control these 
systematic errors, (i) We form a variational basis 
with different junction locations (x, y in Fig. I). 
(ii) We use multihit for the timelike links. (Hi) 
We generalize the multilevel algorithm of @ , orig- 
inally proposed for Polyakov loop correlators, to 
baryonic Wilson loops. This method provides a 
variance reduction exponential in T, which allows 
us to extract the potential from longer loops, with 
crucially improved filtering of excited states. 

3. Results 

A sample of current results based on 160 ana- 
lyzed 1 6 3 x 32 configurations at /3 = 5.8 and 6.0 is 
shown in Fig. 2 (3 quarks in an equilateral trian- 
gle). They are compatible with our earlier mea- 
surements , but the reduced errors now clearly 
show that neither ansatz gives a proper descrip- 
tion of the potential. It approaches the A-ansatz 
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Figure 2. Static potential V qqq vs quark separa- 
tion at f3 = 5.8 and 6.0. Also shown are the A- 
and V-predictions Eqs.(|]) and (||). 

at short distances as expected, but seems to rise 
faster, perhaps as fast as the V-ansatz, at large 
distances. Furthermore, the larger the quark sep- 
aration, the more our variational groundstate fa- 
vors junctions located near the Steiner point. 

To elucidate the asymptotics of the potential, 
we turned to the 3-state Potts model. This toy 
model preserves the center degrees of freedom 
of SU(3) and is thus more likely to agree with 
center-vortex-based predictions of the A-ansatz. 
In this model, we measured the 3-spin correla- 
tion, after adjusting the coupling to match the 
flsu(3) = 6.0 correlation length. High-precision 
cluster Monte Carlo results were obtained for 
multiple 3-spin geometries, in 2c? and 3d. In all 
cases, the 3-spin correlation behaved just like in 
SU(3), falling "in-between" the A- and the Y- 
ansatze. But we could establish that the poten- 
tial was rising asymptotically oc Ly. Large sep- 
arations are required to see this. An example is 
shown in Fig. 3, where the change in action den- 
sity caused by the 3 sources (a) is compared with 



3 




(a) (b) (c) 

Figure 3. Action densities in the 2d Potts model: 
(a) qqq, (b) qq; (c) A prediction, i.e. superposi- 
tion of 3 qq densities. The qq distance is ~ 2.8 fm. 

the superposition (c) of qq action densities (6), as 
predicted by the A-ansatz. A Y string pattern is 
visible, but the strings are "fat" (~ 0.8 fm) com- 
pared to their length, even though the "quark" 
separation is ~ 2.8 fm. 

Remarkably, the approach to Y-asymptotia is 
well described by the same ansatz in the 2d-, 3d- 
Potts model and in SU(3). As shown in Fig. 4, 
the measured value of the 3-spin correlation (or 
Vqqq) falls short of the Y-prediction by an amount 

Vg g q - V q Y qq « c + ci exp(-L A /c2) (3) 

where Co, ci, C2 depend on the qqq triangle geom- 
etry. For the SU{2>) equilateral case, C2 ~ 2 fm. 

4. Conclusions 

Our results show that the baryonic static po- 
tential is neither of the A- nor of the Y-type. It 
approaches the A-ansatz at short distances, but 
rises like the Y-ansatz at large distances. For an 
equilateral qqq arrangement, departure from the 
A-ansatz is not significant until d qq ~ 0.7 fm, so 
that the A-ansatz may be the more relevant one 
for quarks confined inside a hadron. 

The delay in the onset of the Y-behavior is pre- 
sumably caused by fluctuations in the location of 
the junction. In the Potts model, locations xj 
away from the Steiner point are suppressed only if 
the associated length L(xj) of the "glue" strings 
exceeds the minimum Ly by an amount compa- 
rable to the correlation length £. For a junction 
coinciding with one of the "quarks" , for example, 
the inequality L(xj) — Ly > £ becomes satisfied 
only if d qq > 1.7 fm. 

Unfortunately, the Y-behavior at large dis- 
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Figure 4. Amount by which V qqq misses the Y- 
prediction vs perimeter of the quark triangle, for 
the 3d Potts (top) and SU (3) (bottom) cases. 

tances does not help to rule out the center vortex 
picture of confinement. Rather, the failure of the 
A-type prediction of |3| appears to be due to ad- 
ditional hidden assumptions about the indepen- 
dence of certain linking numbers. 
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